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Abstract 

The polarized lepton pair forward-backward asymmetries in B — > p£ + £~ decay 
using a general, model independent form of the effective Hamiltonian is studied. The 
general expression for nine double-polarization forward-backward asymmetries are 
calculated. The study of the forward-backward asymmetries of the doubly-polarized 
lepton pair proves to be very useful tool in looking for new physics beyond the standard 
model. 

PACS numbers: 13.20.He, 12.60.-i, 13.88. +e 



*e-mail: taliev@metu.edu.tr 
^e-mail: savci@metu.edu.tr 



1 Introduction 



Rare B meson decays, induced by flavor changing neutral current (FCNC) b — > s(d)£ + £~ 
transitions provide a promising ground for testing the gauge structure of the Standard 
Model (SM). These decays which are forbidden in the SM at tree level, occur at loop level 
and allow us to check the prediction of the theory at quantum level. Moreover, these decays 
are also quite sensitive to the existence of new physics beyond the SM, since loops with new 
particles can give contribution to these decays. The new physics effects in rare decays can 
appear in two ways; one via modification of the existing Wilson coefficients in the SM, or 
through the introduction of some new operators with new coefficients. Theoretical study of 
the B — > X s £ + £~ decays are relatively more clean compared to their exclusive counterparts, 
since they are not spoiled by nonperturbative long distance effects, while the corresponding 
exclusive channels are easier to measure experimentally. Some of the most important exclu- 
sive FCNC decays are B — > K*^ and B — > (n, p, K, K*)£ + £~ decays. The latter provides po- 
tentially a very rich set of experimental observables, such as, lepton pair forward-backward 
(FB) asymmetry, lepton polarizations, etc. Various kinematical distributions of such pro- 
cesses as B -> K(K*)£+£- [1-3], B -> n(p)£ + £- [4], B s4 -> £+£~ [5] and B s4 -> 1 £ + £~ [6] 
have already been studied. Study of the forward-backward asymmetry, single polarization 
asymmetry, etc., which are experimentally measurable quantities, is an efficient approach 
that have been already been investigated in detail for the B — > K(K*)£ + £~ decay in [1, 7- 
12] in fitting the parameters of the SM and put constraints on new physics. It has been 
pointed out in [13] that some of the single lepton polarization asymmetries might be quite 
small to be observed and might nor provide sufficient number of observables in checking the 
structure of the effective Hamiltonian. In [10] the maximum number of independent ob- 
servables are constructed by taking both lepton polarizations into account simultaneously. 
It is clear that, measurement of many more observables which would be useful in further 
improvement of the parameters of the SM probing new physics beyond the SM. It should 
be noted here that both lepton polarizations in the B — > K*t + t~ and B — > K£ + £~ decays 
are studied in [14] and [15], respectively. As has already been noted, one efficient way of 
establishing new physics effects is studying forward-backward asymmetry in semileptonic 
B — > K*£ + £~ decay, since, Afb vanishes at specific values of the dilepton invariant mass, 
and more essential than that, this zero position of Afb is known to be practically free of 
hadronic uncertainties [12]. The decays of B mesons induced by the b — > d£ + £~ transi- 
tion are promising in looking for CP violation since the CKM factors V tb V t * d , V ub V* d and 
V c bV* d in the SM are all of the same order. For this reason CP violation is much more 
considerable in the decays induced by b — > d transition. So, study of the exclusive decays 
Bd — > (tt, Pi tj)£ + £~ are quite promising for the confirmation of the CP violation and these 
decays have extensively been investigated in the SM [13] and beyond [4]. 

The aim of the present work is studying the polarized forward-backward asymmetry in 
the exclusive B — > p£ + £~ decay using a general form of the effective Hamiltonian, including 
all possible forms of interactions. Here we would like to remind the reader that the influence 
of new Wilson coefficients on various kinematical variables, such as branching ratios, lepton 
pair forward-backward asymmetries and single lepton polarization asymmetries for the 
inclusive B — > X s ( d )£ + £~ decays (see first references in [11,14,17]) and exclusive B — > 
Kt+l~, p£ + £~, ~f£+e-, 7t£ + £~, p£+£~ [1,2,6,9,18,19] and pure leptonic B -> £+£~ decays 



1 



[5,20] have been studied comprehensively. 

The paper is organized as follows. In section 2, using a general form of the effective 
Hamiltonian, we obtain the matrix element in terms of the form factors of the B — > p 
transition. In section 3 we derive the analytical results for the polarized forward-backward 
asymmetry. Last section is devoted to the numerical analysis, discussion and conclusions. 

2 Calculation of double lepton polarizations in B — > 

p£ + £~ decay 

In this section we calculate the double lepton polarizations using a general form of the 
effective Hamiltonian. The B — > p£ + £~ process is governed by b — > d£ + £~ transition at 
quark level. The matrix element for the b — > d£ + £~ can be written in terms of the twelve 
model independent four-Fermi interactions in the following form: 

Heff = ^VuvAcsLdia^Lbhn + CBRdia^Rbhn 
V27T I q q 

+ ClL d L i,b L h-fh + c%k d L i,b L l Rl n R + c RL d Rl ,b R l Ll n L 

+ C RR djfi^bn IbhHr + C LRLR d L b R £ L £ R + C RLLR d R b L £ L £ R 
+ C LRRL d L b R l R l L + C R lrl d R b L £ R £ L + C T da^b la^l 

+ tC TE e^da^b£a a ^ , (1) 

where L and R in (1) are defined as 

T 1-75 n 1+75 



R = 



1 r, ' 
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and Cx are the coefficients of the four-Fermi interactions. The first two coefficients in Eq. 
(1), C SL and C BR , are the nonlocal Fermi interactions, which correspond to -2m s C e 7 ff and 
—2m b Cj^ in the SM, respectively. The following four terms with coefficients C LL , C LR , C RL 
and C RR are the vector type interactions. Two of these vector interactions containing C L °[ 
and C L ot R do already exist in the SM in the form (Cg ff - C w ) and (Cg ff + C w ). Therefore, 



C L °[ and C L ° R can be written as 

^LL — °9 — U 10 + ^LL , 
^LR — ^9 + °10 -T^LR, 

where Cll and Clr describe the contributions of the new physics. The terms with coeffi- 
cients Cl R l R , C R ll R , Cl RR l and C R l R l describe the scalar type interactions. The remain- 
ing last two terms lead by the coefficients Ct and Cte, obviously, describe the tensor type 
interactions. 

It should be noted here that, in further analysis we will assume that all new Wilson 
coefficients are real, as is the case in the SM, while only Cg^ contains imaginary part and 
it is parametrized in the following form 

C e / f = 6 + Aufc , (2) 
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where 



x V ubVud 

u v tb v: d ' 

and 

6 = 4.128 + 0.138cu(s) + g(m c , s)C (m b ) - ^g(m d , s)(C 3 + C 4 ) 

- ±g(m b , s)(4C 3 + 4C 4 + 3C 5 + C 6 ) + ^(3C 3 + C 4 + 3C 5 + C 6 ) , 
6 = [g(rh c , s) - g(rn u , s)](3Ci + C 2 ) , (3) 

where rh q = m^/m^, s = q 2 , Co(/i) = 3Ci + C2 + 3C3 + C4 + 3C5 + 6*6, and 

= --7T 2 - -Li 2 (s) - - ln(s) ln(l - s) - f + 45 . ln(l - s) 
v ; 9 3 V ; 3 V ; V ; 3(1 + 2s) V ; 

2s(l + s)(l — 2s) 5 + 9s -6s 2 

3(1 - s) 2 (l + 2s) { } 3(1 - s)(l + 2s) ' 1 J 

represents the 0(a s ) correction coming from one gluon exchange in the matrix element 
of the operator 0$ [21], while the function g(m q , s) represents one-loop corrections to the 
four-quark operators Oi~O e [22], whose form is 

^ 8, , A > 8 4 2, n 
g(m q , s) = -- \n(m q ) + — + -y q - - (2 + j/,) 



)/|i-y,l{e(i-y ff ) m(y 



1 + ^/1=^ . ' 

1 m 



+ 6{y q - 1) arctan y ^-^j j j - ( 5 ) 



where = 4m 2 /s. 

In addition to the short distance contributions, B — > X d £ + £~ decay also receives long 
distance contributions, which have their origin in the real mm, dd and cc intermediate states, 
i.e., p, u and J/ip family. There are four different approaches in taking long distance 
contributions into consideration: a) HQET based approach [23], b) AMM approach [24], 
c) LSW approach [25], and d) KS approach [26]. In the present work we choose the AMM 
approach, in which these resonance contributions are parametrized using the Breit-Wigner 
form for the resonant states. The effective coefficient including the p, u and J/ip 

resonances are defined as 

C e 9 ff = C 9 (p) + Y res (s) , (6) 

where 



Vi=p,ui Vi ^ *i v% ) 
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The phenomenological factor Ki has the universal value for the inclusive B — > X s ^)£ + £~ 
decay Ki ~ 2.3 [27], which we use in our calculations. 

The exclusive B — > p£ + £~ decay is obtained from the matrix elements of the quark 
operators in Eq. (1) over meson states, which can be parametrized in terms of the form 
factors. Obviously, the following matrix elements 



B 



dia pu q v (l±^)b 



B 



d(l± 75 )6 
da^b B 



B 



are needed in obtaining the decay amplitude of the B 
elements are defined as follows: 



p£ + £ decay. These matrix 



Pip p , e) |^7m(! ± 7s)&| B(p B )) 
(7 2V(^ 



trig + m p 



±i£*(m B + m p )A 1 (q ) 



(8) 



=Fi(Pfl+P,)p(e*g) M f ] =F [^(<7 2 ) - A (g 2 ) 

p(p p ,e) dia^il ±75) ft £(p B 

2 



4e^ ACT 5*>yr 1 (g 2 )±2z e ;(m| - m 2 ) - (p B + P^U^) T 2 (q 2 ) 



(9) 



±2i(e*g) 



Q 



m% — m 2 



P(P P , e) 



5(pb)> = 



ie^ A J - 2T 1 (< 7 V A (p B + + -^(m B - m 2 ) T^q 2 ) - T 2 (q 2 ) 



e™q u 



(10) 



T^q 2 ) - T 2 {q 2 ) - 



m, 



rat 



r% 2 ) 



where q = Pb~P p is the momentum transfer and £ is the polarization vector of p meson. In 
order to ensure finiteness of (8) and (10) at q 2 = 0, we assume that A 3 (q 2 = 0) = A (q 2 = 0) 



and Ti(q 2 = 0) = T 2 (q 2 = 0). The matrix element (^p d(l ± 75)6 
contracting both sides of Eq. (8) with q^ and using equation o 
mass of the d quark we get 



Bj can be calculated by 
motion. Neglecting the 



p(p p ,e) d(l±7 5 )6 B(p B ) 



1 

m b 



T'2im p ( y e*q)A (q 2 



(11) 



In deriving Eq. (11) we have used the relationship 

2m p A 3 (q 2 ) = (m B + m p )A 1 (q 2 ) - (m B - m p )A 2 (q 2 ) , 
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which follows from the equations of motion. 

Using the definition of the form factors, as given above, the amplitude of the B — > pi + t 
decay can be written as 

M{B^ P tn = ^ f v tb v; d 



x I - I - 2A 1 e^ Xa E*» P y ~ iB^l + iB 2 (e*q)(p B + Pp ), + iB 3 (e*q)q^ 



+£-f(l + [ - 2C 1 6^ a e* l/ P y - iD 1 e; + iD 2 (e*q)(p B + Pp )^ + tD 3 (e*q)q p 
+1(1 - lh )l[iBt{e*q)] + 1(1 + -y 5 )e[iB 5 (e*q) 

+Ua^t(iC T e^) [ - 2T 1 e*\p B + Pp Y + B 6 e* x q° - B 7 (e*q)p p Y 



+ 16C TE Ia p J - 2T 1 e^(p B + Vp T + B % e**<? - B 7 (e*q) Pp »q v , 



(12) 



where 



B 1 
B 2 

B 3 

Ci 
D 1 
D 2 
D 3 
B, 

B 5 
B 6 
B 7 



V T 
(C% + Crl)^^- - 2(C BR + C SL )^ , 

(Ci° L - C RL )(m B + m p )A x - 2(C BR - C SL )(m 2 B - m 2 p )% , 



/-(tot /~i 1 

CLL ~^ L A 2 -2(C BR -C SL )1 
m B + m p q l 



T 2 + 



— 2 9^3 



^3-^0 



2(C^ - C RL )m p 

q z 

s~itot s~i 
^LR > 
itot 



2~3 

+ 2(C BR — Csl)—^ , 



tot 
LL 



c 



RR) 



Ai(C 

r> / s-itot . ritot r~i . /m \ 

&l{ U LL ~^ ^LR > ^RL — >■ ^i?i?J , 

d ( s-itot . rttot /~i . si \ 

HiyyLL ~* ^LR > ^RL —> L'RR) , 

73 / sitot . sitot si . si \ 

a A^LL ~^ ^LR > U RL — >• ^i?i?J > 

—2(Clrrl — C R l R l) — ~A) ; 
—2(C B rlr — Crllr) — -Aq , 
2(m 2 B -m 2 p ) 1 - 2 , 



m n — m; 



(13) 



From this expression of the decay amplitude, for the differential decay width we get the 
following result: 



§(5 - pt£-) = \V tb V;/ A 1/2 (l, f, s)vA(s) , (14) 
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with 



3f p s 



m% Re[ - 6m B m e s\(B 1 - D^Bl - B* 5 ) 



- l2m 2 B rh 2 e s\{B 4 B* + (B 3 - D 2 - D Z )B{ - (B 2 + £? 3 - D 3 )Dl} 
+ Qm%m e s(l - r p )X(B 2 - D 2 )(B* A - B* 5 ) 

+ 12m B m 2 s(l - r p )X(B 2 - D 2 )(B* 3 - D* 3 ) 
+ 6m 3 B m e Xs 2 (B 4 - B b )(B* z - £)*) 
+ 48m^f p s{3fiiDt + 2m^AA 1 C 1 *} 
+ 48m|m^A 2 ( J B 2 + D 2 )B* 7 C* TE 

- lQm B f p s(m 2 - s)A{ \A X \ 2 + IC^ 2 } 

- m|s(2m 2 - s)A{ |5 4 | 2 + \B 5 \ 2 } 

- 48m|m £ s(l - r p - s)A{(Bi + D X )B*C* TE + 2(5 2 + £) 2 )S*CJ B } 

- 6m^m 2 sA{2(2 + 2f p - s)5 2 D 2 * - s |(B 3 - L> 3 )| 2 } 
+ 96m B m £ s(A + 12f p s)(£i + D^)BIC* TE 

+ 8m|s 2 p |C T | 2 + 4(3 - 2v 2 ) \C TE \ 2 }{4(A + 12f p s) \B & \ 2 

- 4m|A(l -r p - s)B e B* 7 + m^A 2 |5 7 | 2 } 

- 4m|A{m 2 (2 - 2f p + s) + s(l - r p - s)}( J B lJ B 2 * + 
+ s{6f p s(3 + v 2 ) + A(3 - v 2 )}{ Ifiil 2 + |A| 2 } 

- 2m^A{m 2 [A - 3(1 - f p f] - As}{ \B 2 \ 2 + \D 2 \ 2 } 
+ 128m|{4m 2 [20f p A - 12f p (l - f p f - As] 

+ s[4f p A + 12f p (l - r p f + AS]} |C T | 2 |Tx| 2 

+ 512m|{s[4f p A + 12f p (l - f p f + As] 

+ 8m 2 [12f p (l - f p f + A(s - 8f p )]} \C TE \ 2 \T,\ 2 

- 64m|s 2 p \C T \ 2 + 4(3 - 2v 2 ) \C TE \ 2 }{2[A + I2f p {\ - r p )]B 6 T* 

- m 2 B X(l + 3f p -s)B 7 T*] 

+ 768m 3 B m e r p sX(A 1 + d)C^T* 

- 192m B m^s[A + 12f p (l - f p )](Bi + D 1 )C^ E T* 

+ 192m|m^sA(l + 3f p - s)A(5 2 + D 2 )C; e T*] , (15) 

where s = (f /m 2 B , f p = m 2 /m 2 B and A(a, 6, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 2bc, rh^ = mi/m B , 

v = \Jl — Am 2 1 s is the final lepton velocity. 

Using the matrix element for the B — > p£ + £~ decay, our next problem is to calculate 
the polarized FB asymmetries. For this purpose, we define the following orthogonal unit 
vectors sf M in the rest frame of i^, where % = L, N or T correspond to longitudinal, normal, 
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transversal polarization directions, respectively (see also [1,8,10,15]), 



Sr = 



Sj-< 



S L 



b N — 



Sj-< 



0,e N 
0,ef 

0,e + 
0,e + 



0,^ 

v \P- 



\Pk x p 
(0, x e£" 



) 



°'P1 



0, x e£) 



(16) 



where p T and px are the three-momenta of the leptons £ T and p meson in the center of 
mass frame (CM) of l~ £ + system, respectively. Transformation of unit vectors from the rest 
frame of the leptons to CM frame of leptons can be done by the Lorentz boost. Boosting 
of the longitudinal unit vectors s^ yields 



L J CM 



\Pt\ E ePT 



m e m E \p T \ 



;i7) 



where p + = —p~, Eg and me are the energy and mass of leptons in the CM frame, respec- 
tively. The remaining two unit vectors s^, are unchanged under Lorentz boost. 

The definition of the unpolarized and normalized differential forward-backward asym- 
metry is (see for example [28]) 



f 1 ^ r dz d * L dz 

Jo dsdz J-i dsdz 

r 1 d 2 T ^ | ro 
Jo dsdz J-i 



d 2 T 



Afb — 



o d 2 T , ' 

dz 



;is) 



dsdz 



where z = cos 9 is the angle between B meson and l~ in the center mass frame of leptons. 
When the spins of both leptons are taken into account, the Afb will be a function of the 
spins of the final leptons and it is defined as 



Ap B (s) 



' dT{s) Y 
d 2 T(s, s- 



dz — 



dz 



d 2 T(s,s =i,s + = j) d 2 T(s,s =i,s + = —j) 



ls + = j) d 2 r(s,s- 



dsdz 

— * 

—i, s' 



dsdz 



-1) 



dsdz 



dsdz 



= A F b{s =i,s + = j) 
+ A FB (s- 



A 



fb(s 



-j) - A F b(s' 



■ ->+ 
-t, s 



-3 



(19) 



Using these definitions for the double polarized FB asymmetries, we get the following 
results: 
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A F L B = —rn%\f\v Re[ - m 3 B m e \U(B 1 - D X )B*C* T - {B A + B b )(B* 2 + D*)\ 

+ 4m B m e \{(l - r p )(B 2 - D 2 )B* 7 C* T + s(B 3 - D 3 )B* 7 C* T ) 

- m e (l -r p - s){Bl(B 4 + B 5 - SB 6 C t ) + D*(B 4 + B 5 + 8B 6 C T )} 
+ BmefpsiA^* - dD*) + 128m|m*7>s(Ai - d)B^ E 

+ 2m|sA{(5 4 - B 5 )B* 7 C^ + 2(5 4 + 5 5 )5 7 *C^} 

- 8m|m,(l - r>)(l - f p - s)(S 2 - L> 2 )5 6 *C; 

- Am B {\ -f p - s)s{(B 4 - B 5 )B* 6 C^ + 2{B 4 + B 5 )B* 6 C^ E 

+ 2m B me(B 3 ~ D 3 )B* 6 C* T } - 256m B m e r p (l - r p ){A x - C^T*^ 

- 16m £ (l - 5f p - s)(B 1 - D 1 )T*C^ r 

+ 16m|m,(l - f p )(l + 3f p - s)(B 2 - D 2 )T*C* T 

+ 8m B (l + 3f p - s)s{2(5 4 + B 5 )T*C^ E + (5 4 - 5 5 )T*C* 

+ 2m B m,(fl 3 - £> 3 )7?C?.}] , 
8 

*4fb = mlv^A^Im - m^B-yDl + m B XB 2 D 2 ) + Am^mfi p \TsA x Cl 

oT pS l\ 

- 2m B s{B 6 (C T - 2C TE )B{ + B 6 (C T + 2C TB )^} 

- m%s\{B 7 (C T - 2C TE )B* 2 + B 7 (C T + 2C TE )D* 2 } 

- 16m|m^(4 \B 6 \ 2 + m B \ \B 7 \ 2 )C T C^ E 
+ m|m € (l - r p - s)(B 1 D* 2 + B 2 D{) 

+ m|s(l - r p - s){(S*S 7 + 2B*B 6 )(C T - 2C TE ) 

+ {D\B 7 + 2D*B 6 )(C T + 2C TE )} 

- 64m 2 B m e s{ - m|(l - r p - s)Re[B 6 B;] + 4 |Ti| 2 - 4Re[B 6 T*} 
+ 2m|(l + 3f p - s)Re[B 7 T*]}C T C* TE 

+ 16m|f p s|(Ai - CJC^T* - 2{A 1 + Ci)C^T*} 
+ Am B s{Bl{C T - 2C TE )T X + D\{C T + 2C r s)Ti} 

- 4m|s(l + 3f p - s){s 2 *(C T - 2(7^)71 + L> 2 *(C T + 2C TB )Ti}" , 

^| = mlv^A^Im^ - rh^BxDl] + m%XB 2 D 2 ) + Am B rh(.r p sA 1 C* 1 

oT pS ZA 

+ 2m B s{ J B 6 (C T + 2C TE )5* + £ 6 (C T - 2CW)^} 
+ m|sA{s 7 (C7 T + 2C TE )B* + B 7 (C T - 2C T£ )D;} 
+ 16m|m € s(4 l^] 2 + m^A \B 7 \ 2 )C T C* TE 



+ m%m t (l -f p - s^D; + B 2 Dl) 

- m%s(l -f p - s){(BlB 7 + 2B* 2 B 6 )(C T + 2C TE ) 

+ {D\B 7 + 2D*B 6 )(C T - 2C TE )} 

+ 6Am 2 B m e s{ - m|(l - f p - s)Re[5 6J B;] + 4 |Tx| 2 - 4Re[S 6 T*] 
+ 2m|(l + 3f p - s)Re[B 7 T*]}C T G 



A 



LT 
FB 



+ \Qm\f p s\{Ai - C^T* + 2{A l + C^C^T*} 
- Am B s{Bl{C T + 2Cte)Tx + D{{C T - 2C TE )T 1 } 
+ 4m|s(l + 3f p - s){B* 2 (C T + 2C TE )T 1 + L> 2 *(C T - 2C T b)7i} 

,|v^A Re[ - %{ | Si + £>i| 2 + m^A |B 2 + £> 2 | 2 } 



(22) 



3f p sA 



+ Am B fh e f p s{ \A X + C 1 \ 2 } 

- 64m 2 B m e s \C TE \ 2 {4 |^ 6 | 2 + m B \ \B 7 \ 2 - 4m|(l - f p - s)S 6J B 7 *} 
+ 2m|m<(l -f p - s)(B 1 + A)^* + D* 2 ) 

+ 2m B {\ -f p - s){4m 2 (2B* 2 B 6 + B{B 7 ){C T + 2C TE ) 

- s(2B*B 6 + B{B 7 ){C T - 2C TE ) 



- Am B {Am 2 [BlB & (C T + 2C TE ) - B & D\(C T - 2C TE ) 

- s[B{B & {C T - 2C TE ) - B 6 Dl(C T + 2C TE )] } 

- 2m B \[Am][B* 2 B 7 (C T + 2C TE ) - B 7 D* 2 (C T - 2C TB ) 

- s[b*B 7 (C t - 2C TS ) - B 7 D* 2 (C T + 2C TS )] } 

- 2m|(l - f p - s){4m 2 (2S 6J D 2 * + B 7 £>j;)(C T - 2C TE ) 

- s(2B 6 D* 2 + B 7 D\){C T + 2C TE )} 

2B 6 T* - m|(l + 3f p - s)S 7 7t 

V| 2 + (4f p -s)|C TE | 2 ]} 

c 1 )c*r* + 2(A 1 -c 1 )c; E r; 



+ 256m 2 3 m£|2s |Ct£ 
+ 4|T 1 | 2 [r p |G 
+ 32m 3 B f p {Am 2 [(A 1 



+ s[a;(C t - 2C TE )T X + cr(C T + 2C TE )T^ } 

+ 8m B {4m 2 (C T + 2C TE ) - s(C T - 2C TE )}{ J B* - m|(l + 3f p - ^B*}^ 
- 8m B {4m 2 (C T - 2C TE ) - s(C T + 2C te )}{D* 1 - m|(l + 3f p - s)^*}^] 



(23) 



A TL 



3f p sA 



m 2 B V~s\ Re[m e { \B 1 + D^ 2 + m%\ \B 2 + D 2 | 2 } 



- Am%m l f p {\A l + C l \ 2 } 
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+ QAm 2 B m e s \C TE \ 2 {4 \B 6 \ 2 + m B \ \B 7 \ 2 - 4m|(l -r p - s)B G B* 7 ) 

- 2m|m/(l - r> - + DO(S* + D* 2 ) 

+ 2m|(l - r p - s){Am 2 (2B* 2 B & + B\B 7 ){C T - 2C TE ) 

- s(2B*B 6 + B*B 7 )(C T + 2C TE ) 



~ 4m B {4m 2 [b\B & {C t - 2C TE ) - B % D\{C T + 2C TE ) 

- s[B\B % {C T + 2C TE ) - B 6 D\{C T - 2C TE )] } 

- 2m%\{Am 2 l [B* 2 B 7 (C T - 2C TE ) - B 7 D* 2 (C T + 2C TE ) 

- s[b* 2 B 7 (C t + 2C TE ) - B 7 D* 2 (C T - 2C TE )] } 

- 2m%(l -f p - s){4m 2 (2B 6 D* 2 + B 7 D{){C T + 2C TE ) 

- s(2B 6 D* 2 + B 7 D\){C T - 2C TE )} 

- 256m|m £ {2s \C TE \ 2 [2B 6 T* - m|(l + 3f p - s)5 7 T* 
+ 4|T 1 | 2 [r p |C T | 2 + (4f p -s) |C TS | 2 ]} 

- 32m 3 B r p {4m 2 [(A 1 + d)C*T* - 2(A X - C^C^T* 
+ s 



A*(C T + 2C TE )T X + cr(C T - 2C TE )T^ } 
- 8m B {4m 2 (CT + 2C TE ) - s(C T - 2C te )}{d{ - m|(l + 3f p - s)^*}^ 
+ 8m B {4m 2 (C T - 2C T£ ) - s{C T + 2C TE )}{ J B* - m|(l + 3f p - s)5 2 *} ,7\ 



A™ 



v / AIm[m|m^A{( J B 4 - S 5 )(S^ + D* 2 ) + %B 7 C TE (B{ - D{) 



+ 8m 2 B sB;CT E (B 3 - A 
- 2m B m 2 s\(B 2 + D 2 )(B* 3 - £>*) 

+ 4m^m,(l - f p )A{2m i? s J B 7 *C* s ( J B 2 - L> 2 ) + m^D*} 
+ 2m 2 B m 2 e s(l + 3f p - s)(5i J B 2 * — DiD 2 ) 



E*(5 4 - E 5 + 165 6 C TB ) 



- L>?(5 4 - B 5 - 16B G C TE ) + 2m B m e (Bi + D 1 )(B* 3 - D* 3 ) 



+ 4 



rh e B 1 D* 1 + Am 3 B s 2 B 6 C TE (B* 3 - D* 3 )] } 



- lGm 3 B m e s(l - f p ){\ -f p - s){B 2 - D 2 )B* 6 C* TE 
+ 2m 2 B m 2 [\ + (1 - f p )(l - r> - s)](SJD 2 + B^Di) 
+ 32m|m/s(l - r p )(l + 3r> - s)(B 2 - D 2 )C* TE T* 

- 8m B s(l + 3f p - s){4m £ (5! - DJC^T* - 2m B s(B A - B 5 )C^ E T* 

- Am 2 B m e s(B 3 - D Z )C* TE T* + m B sv 2 (B 4 + 5 5 )C£T*} 
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- 4m|s 2 (l - r p - s){2{B A - B 5 )B* 6 C TE - v 2 (B 4 + B 5 )B* 6 C^} 

+ 2m B s 2 \{2(B 4 - B 5 )B;C^ E - v 2 {B 4 + B 5 )B* 7 C^}] , (25) 



A NN _ _ A TT 
M-FB — ^-FB 



r p A 



m 3 B V\v Re[ - m 2 B m e \{A{B 1 - D X )B*C* T + (B 2 + D 2 ){Bl + B* b )) 



+ Am B rh e X{(l - r p ){B 2 - D 2 )B*C* T + s(B 3 - D 3 )B* 7 C* T } 
+ 2m|sA{(5 4 - B b )B*C^ - 2(B 4 + B b )B* 7 C* TE ] 
+ m e (l -f p - s){Bl(B 4 + B 5 + 8B 6 C T ) 
+ Dl(B 4 + B 5 -8B 6 C T )} 

- %m 2 B m t {\ - r p ){\ -f p - s)(B 2 - D 2 )B;C* T 

- Am B s(l -r p - s){(B 4 - B 5 )B* 6 C^ - 2(B 4 + B 5 )B* 6 C* E 

+ 2m B rh e (B 3 - D 3 )B* 6 C^} 

+ lGm%m e (l - r p )(l + 3r p - s)(B 2 - D 2 )C*T* 

+ 8m B s(l + 3f p - s){(B 4 - B 5 )C*T* - 2(B 4 + B 5 )C^ E T*} 

- 16m £ (l + 3f p - s)(Si - D^C^T* 

+ lQm 2 B m e s(l + 3f p - s)(B 3 - D 3 )CtT*] . (26) 

In these expressions for A EB , the first index in the superscript describes the polarization 
of lepton and the second index describes that of anti-lepton. 

3 Numerical analysis 

In this section we analyze the effects of the Wilson coefficients on the polarized FB asymme- 
try. The input parameters we use in our numerical calculations are: m p = 0.77 GeV, m T = 
1.77 GeV, m M = 0.106 GeV, m b = 4.8 GeV, m B = 5.26 GeV and T B = 4.22 x 10~ 13 GeV. 
For the values of the Wilson coefficients we use Cf M = -0.313, Cf M = 4.344 and 
Cf M = —4.669. It should be noted that the above-presented value for Cg M corresponds 
only to short distance contributions. In addition to the short distance contributions, it 
receives long distance contributions which result from the conversion of uu, dd and cc to 
the lepton pair. In order to minimize the hadronic uncertainties we will discard the regions 
around low lying resonances p, w, J/ip, if)', ip", by dividing the q 2 region to low and high 
dilepton mass intervals: 

Region I: 1 GeV 2 < q 2 < 8 GeV 2 , 

Region II: 14.5 GeV 2 < q 2 < (m B - m p ) 2 , 

where the contributions of the higher ip resonances do still exist in the second region. For 
the form factors we have used the light cone QCD sum rules results [24,29]. As a result 
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of the analysis carried out in this scheme, the q 2 dependence of the form factors can be 
represented in terms of three parameters as 



F(0) 



1 — dpS + &fs 2 
where the values of parameters -F(O), ap and bp for the B 



p decay are listed in Table 1. 





F(0) 

V"/ 


(Xp 


bp 




0.372 ±0.04 


1.40 


0.437 


A B r p 


0.261 ±0.04 


0.29 


-0.415 


^2 


0.223 ±0.03 


0.93 


-0.092 


yB^p 


0.338 ±0.05 


1.37 


0.315 


T B^p 


0.285 ±0.04 


1.41 


0.361 


rpB-^p 


0.285 ±0.04 


0.28 


-0.500 




0.202 ±0.04 


1.06 


-0.076 



Table 1: B meson decay form factors in a three-parameter fit, where the radiative correc- 
tions to the leading twist contribution and SU(3) breaking effects are taken into account. 



In further numerical analysis, the values of the new Wilson coefficients are needed, and in 
the present analysis we will vary them in the range — |Cio | < |Cj| < |C 10 | . The experimental 
value of the branching ratio of the B — > K*£ + £~ decay [30,31] and the bound on the 
branching ratio of the B — > p + p~ [32] suggest that this is the right order of magnitude for 
the vector and scalar interaction coefficients. It should be noted here that the experimental 
results lead to stronger restrictions on some of the Wilson coefficients, namely —1.5 < Cp < 
1.5, —3.3 < Cte < 2.6, —2 < Cll, Crl < 2.3, while the remaining coefficients vary in the 
range —4 < Cx < 4. 

As is obvious from the explicit expressions of the forward-backward asymmetries, they 
depend both on q 2 and the new Wilson coefficients Cx- As a result of this, it might be 
difficult to study the dependence of the polarized forward-backward asymmetries A FB on 
these parameters simultaneously. Therefore, it is necessary to eliminate the dependence 
of A l p B on one of the parameters. We eliminate the dependence of the polarized A FB 
on q 2 by performing integration over q 2 in the kinematically allowed region, so that the 
polarized forward-backward asymmetry is said to be averaged. The averaged polarized 
forward-backward asymmetry is defined as 



A ij 



FB 



(I 2 

■irtia-.r 



2 

min 



A% B f~ 2 dq 2 
dq 2 



I 



2 

I max 



<m 

dq~ 2 



^dq 2 
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In Figs. (1) and (2), we present the dependence of \ApBj 011 Cx for the B — > pp + p~ 
decay, in the Regions I and II, respectively. The common intersection point of all curves 
corresponds to the SM case. We observe from this figure that, (^Ap B ^ has practically 
symmetric behavior in regard to its dependence on Cp and Cpp with respect to zero position. 
We see from Fig. (1) that, (Ap 1 ^ for the B — > pp + p~ decay is very strongly dependent 
on Cll, Clr, Cp, Cpp and on Crl when Crl > 1. There are certain regions of Cx where 
the magnitude of (Ap 1 ^ is, more or less, two times larger than compared to its value in 
the SM. This fact is a direct indication of the confirmation of new physics beyond the SM 
which can be attributed to the existence of new vector type interaction. (Ap L B ^ behaves 
in the same way in region II as it does in Region I, except for the scalar interactions and 
vector interaction with coefficient Crr (see Fig. (2)), and being quite sensitive to the rest 
of the remaining new Wilson coefficients. It is interesting to observe that the sign of (Ap 1 ^ 
is negative (positive) in Region I (Region II) for all values of Cx- 

Figs. (3) and (4) depict the dependence of {Ap T ^ on Cx for the B — > pp + p~ decay, in 
Regions I and II, respectively. We observe from Fig. (3) that, except scalar interactions, 
{Ap 1 ^! is quite sensitive to the existence of the remaining ones. More important than that 

is the presence of regions of the new Wilson coefficients where {Ap 7 ^ changes its sign, while 

in the SM case the sign of the (A F T ^} is never switched. So, study of the magnitude and 

sign of \A F T B j can serve as a good test for looking new physics beyond the SM. In region 

II, (AJpB*} is strongly dependent only on tensor interactions (see Fig. (4)). 

The dependence of (^Ap B ^ on Cx for the B — > pp + p~ decay, is given in Fig. (5) in 

Region I, and Fig. (6) in Region II, respectively. We observe from these figures that (^A FB ^ 
exhibits strong dependence only on tensor interactions, and especially there is a region of 
C T e where {Aps^ exceeds the SM prediction more than one order of magnitude. Moreover, 

when Cp and Cpp is negative (positive) the sign of (A^/p^j is positive (negative). Hence, 
determination of its magnitude and sign is an unambiguous confirmation of the existence 
of tensor interaction. Similarly, Fig. (6) depicts strong dependence of (Ap^p^j on tensor 

interactions. When Cpp is negative the sign of (^Ap^j is positive, and when Cpp is positive 
it is negative. Analogous behavior is observed for Cp, except for quite a narrow region. 

All remaining forward-backward asymmetries for the B — > pp + p~ decay are numerically 
very small and for this reason we do not present them. 

The study of the dependence of forward-backward asymmetry for the B — > pr + T~ 
decay gives richer information. In Figs. (7), (8), (9) and (10) we present the dependence 
of (A L F L B ), (A%), (A t fb ), {A N F T B ) = (A T F N B ) and (A™) = - (A T F T B ) on new Wilson 
coefficients, respectively. We see from Fig. (7) that when C x is negative, (A LL ^ > (Asm) 
only for the scalar interaction C RLLR . Also, when C x is positive, (A LL ^ > (Asm) f° r 

the coefficients Crl, Clrlr. Fig. (8) depicts that (A F T B ^ is strongly dependent on tensor 
interactions, as well as on vector interactions with the coefficients Cpp and Crr when 
they get positive values. We observe from Fig. (9) that, (^A FB ^ = {A^q^ both exhibit 
strong dependence on all new Wilson coefficients. As can easily be observed from Fig. (10), 
(A fb ^ = — (^ApT^j both are very sensitive to the presence of tensor and scalar interactions. 
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It follows from these results that few of the polarized forward-backward asymmetries 
show considerable departure from the SM predictions and these ones are strongly dependent 
on different types of interactions. Hence, the study of these quantities can play crucial role 
in establishing new physics beyond the SM. 

At the end of this section, we would like to discuss the following problem. It is clear 
that the existence of new physics can more easily be checked through branching ratio 
measurements. In this connection there follows the question: could there be a situation in 
which the branching ratio coincides with that of the SM result, while polarized forward- 
backward asymmetry does not? In order to answer this question we study the correlation 



between the (Jip B j and the branching ratio B. We can briefly summarize the results of 
our numerical analysis as follows: As far as B — > pfi + fi~ decay is concerned, except for a 
very narrow region of Crr, such a region is absent for all new Wilson coefficients for the 
asymmetries (Ap 1 ^), \A F b) and \A FB ). 



The B — > pr + r decay is more informative for this aim. In Figs. (11) and (12) we 
present the dependence of (Ap^j and (Ap T B ^j on the branching ratio. It follows from these 
figures that, there indeed exists certain regions of Cx for which the polarized forward- 
backward asymmetry differs from the SM prediction, while the branching ratio coincides 
with that of the SM result. We also note that, such a region exists for the asymmetries 
\Apb) = — \A FB ) as well, only for the tensor interaction. 



In conclusion, in this work we investigate the forward-backward asymmetries when both 
leptons are polarized, using a general, model independent form of the effective Hamiltonian. 
We see that the study of the zero position of (A fb ^j can give unambiguous conformation 
of the new physics beyond the SM, since when new physics effects are taken into account, 
the results are shifted with respect to their zero positions in the SM. We also find out that 
the polarized A F b is quite sensitive to the existence of some of the new Wilson coefficients. 
We see that there exist certain regions of some of the new Wilson coefficients for which, 
only study of the polarized forward-backward asymmetry gives invaluable information in 
establishing new physics beyond the SM. 
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Figure captions 



Fig. (1) The dependence of the averaged forward-backward double-lepton polarization 
asymmetry (A F on the new Wilson coefficients C x , for the B — > pp + decay, in Region 
I. 

Fig. (2) The same as in Fig. (1), but in Region II. 

Fig. (3) The same as in Fig. (1), but for the averaged forward-backward double-lepton 
polarization asymmetry \Ap L ^). 



Fig. (4) The same as in Fig. (3), but in Region II. 

Fig. (5) The same as in Fig. (1), but for the averaged forward-backward double-lepton 
polarization asymmetry (Aps). 



Fig. (6) The same as in Fig. (5), but in Region II. 

Fig. (7) The dependence of the averaged forward-backward double-lepton polarization 
asymmetry (AJpfej on the new Wilson coefficients C x , for the B — > pr + r~ decay, in Region 
II. 

Fig. (8) The same as in Fig. (7), but for the averaged forward-backward double-lepton 
polarization asymmetry \A F T B ). 




Fig. (9) The same as in Fig. (7), but for the averaged forward-backward double-lepton 
polarization asymmetry (Ap^) = {A^). 



Fig. (10) The same as in Fig. (7), but for the averaged forward-backward double-lepton 
polarization asymmetry (Ap^) = — \Ap^ ). 



Fig. (11) Parametric plot of the correlation between the averaged forward-backward 
double-lepton polarization asymmetry (Ap 1 ^ and the branching ratio for the B — > pr + T~ 
decay, in Region II. 

Fig. (12) Parametric plot of the correlation between the averaged forward-backward 
double-lepton polarization asymmetry (Ap 7 ^ and the branching ratio for the B — > pr + r~ 
decay, in Region II. 
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